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ABSTRACT 

A  two-dimensional  elliptic  computer  model  was  adapted  to 
the  solution  of  the  flow  field  in  a  turbojet  test  cell 
section  which  contained  the  engine  exhaust  duct  and 
augmenter  tube.   Various  engine  power  settings,  augmentation 
ratios,  and  augmenter  geometries  were  investigated. 

Flow  visualization  and  static  pressure  recovery  were 
primary  objectives  subject  to  specified  engine  power  settings 
and  augmentation  ratios.   Capabilities  and  limitations  of 
the  model  are  discussed. 
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I.   INTRODUCTION 

A.   BACKGROUND 

Jet  engines  undergoing  maintenance  are  tested  in  jet 
engine  test  cells  prior  to  their  operational  use.   A  typical 
test  cell  is  shown  in  Figure  1  and  can  be  viewed  as  a  "U" 
shaped  concrete  structure  of  rectangular  cross-section.   The 
cell  is  designed  to  provide  clean  distortion  free  airflow 
to  the  engine  inlet/  which  is  necessary  for  accurate 
performance  evaluation.   In  addition,  the  cell  attenuates 
noise  and  provides  safety  for  test  equipment  and  personnel. 

Engine  exhaust  is  directed  into  an  augmenter  tube,  then 
deflected  vertically  to  be  expelled  into  the  atmosphere. 
The  augmenter  tube  acts  in  combination  with  engine  exhaust 
momentum  to  form  a  jet  pump.   In  reality,  the  engine  exhaust 
jet  entrains  secondary  air  (air  not  ingested  by  the  engine) 
which  is  then  carried  through  the  augmenter  tube  by  exhaust 
jet  momentum.   Frictional  losses  associated  with  the  mass 
of  air  moving  through  the  test  cell  are  thus  overcome  by  the 
jet  pump's  energy  addition. 

The  quantity  of  secondary  air  entrained  is  crucial  to 
proper  engine  testing.   Too  little  entrainment  allows  exhaust 
gas  recirculation  to  the  engine  inlet.   Too  much  entrainment 
causes  excessive  static  pressure  gradients  to  exist  between 
engine  inlet  and  exhaust  planes.   Static  pressure  should  be 
nearly  identical  in  both  planes  for  proper  testing  but  a 
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small  gradient  must  exist  to  prevent  exhaust  gas  recircula- 
tion.  Sufficient  secondary  air  should  be  entrained  to  cool 
the  engine  exhaust  to  acceptable  levels.   Also,  exhaust  flow 
with  very  high  kinetic  energy  can  damage  the  test  cell; 
secondary  air  entrainment  reduces  this  energy. 

Augmentation  ratio  is  defined  as  the  ratio  of  entrained 
secondary  air  to  the  jet  exhaust  mass  flow  (primary  flow) . 
In  practice,  engine/augmenter  tube  separation  is  adjusted  to 
generate  sufficient  secondary  flow  to  prevent  exhaust  gas 
recirculation  at  all  power  settings.   Occasionally,  the 
augmenter  inlet  diameter  is  varied  to  further  control  the 
entrainment  process. 

Both  analytical  and  empirical  methods  have  been  applied 
to  the  entrainment  mixing  problem.   Empirical  equations  are 
usually  applied  to  particular  test  cells,  but  operator 
experience  coupled  with  trial  and  error  are  required  to 
create  the  desired  secondary  flow. 

Larger  and  more  expensive  test  cells  are  needed  in  the 
future  to  replace  obsolete  cells  and  to  accommodate  new 
engines  with  greatly  increased  mass  flow  rates.   Present 
design  tools  do  not  adequately  describe  the  flow  field  in 
a  test  cell  and  new  design  concepts  cannot  be  evaluated  until 
the  cells  are  constructed.   An  accurate  flow  model  would 
provide  a  needed  design  tool  and  would  help  prevent  costly 
design  errors. 
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B.   TEST  CELL  POLLUTION 

Test  cell  installations  are  becoming  increasingly  sub- 
jected to  localized  state  and  community  pollution  control 
board  regulations.   The  occasional  prominent  exhaust  plume 
rising  from  test  cells  is  a  continual  source  of  local  com- 
plaint.  Even  when  the  plume  is  not  visible,  the  cell  emits 
pollutants  similar  to  those  of  aircraft  operating  in  the 
area. 

Fortunately,  the  stationary  cell  is  amenable  to  exhaust 
treatment  which  has  been  shown  to  be  effective  in  reducing 
particulate  emissions  to  acceptable  levels.   The  cost  asso- 
ciated with  the  construction  and  operation  of  the  necessary 
equipment  is  high  and  is  a  function  of  total  mass  flow 
(primary  plus  secondary)  through  the  cell.   Thus,  it  is 
necessary  to  minimize  total  mass  flow  to  minimize  abatement 
costs. 

Carbon  particulates  and  unburned  hydrocarbons  are  of 
immediate  concern  in  the  abatement  of  emissions  from  jet 
test  cells.   Also  of  concern  are  nitrogen  oxides  (NO  ) , 
carbon  monoxide  (CO)  and  sulfur  dioxide.   The  amount  of  each 
pollutant  emitted  is  a  function  of  engine  type  and  mechanical 
condition.   Power  setting  alters  emissions  with  CO  and  unburned 
hydrocarbons  predominate  at  idle.   Particulates  and  nitrogen 
oxides  are  the  main  pollutants  at  high  power  settings.   After- 
burner operation  provides  very  high  external  exhaust  tempera- 
tures and  chemical  reactions  continue  as  the  flow  passes  into 
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the  augmenter  tube.   The  secondary  air  mixing  is  believed 
to  alter  chemical  reaction  rates  by  effectively  quenching 
portions  of  the  exhaust  jet.   Present  flow  models  do  not 
describe  the  process  completely  and  it  has  not  been  possi- 
ble until  recently  to  accurately  measure  mass  emission  fac- 
tors in  the  afterburner  core  due  to  the  high  temperatures 
and  the  chemical  kinetic  nature  of  the  process. 

To  fully  optimize  the  pollution  abatement  program  it  is 
necessary  to  understand  how  exhaust  pollutants  are  formed, 
in  what  quantities  and  at  what  rates.   The  effect  of  engine 
power  setting  and  secondary  flow  entrainment  on  emitted 
pollutants  must  also  be  analyzed.   Pollutant  measurement 
would  be  simpler  and  less  costly  if  accurate  species  and 
velocity  distribution  predictions  could  be  made  at  specified 
points  in  the  test  cell  exhaust  system.   To  accomplish  these 
tasks,  more  complete  knowledge  of  the  air  flow  through  test 
cells  must  be  acquired. 

C.   SCOPE  OF  STUDY 

Most  analysis  for  predicting  dilution  air  is  very  empiri- 
cal, usually  requiring  an  empirical  jet  spread  equation. 
Recent  advances  have  been  made  toward  solving  the  recircu- 
lating, turbulent  fluid  flow  problem.   One  technique  is  due 
to  Spalding  et   al.  [Refs.  1  and  2].   The  governing  two- 
dimensional  elliptic  partial  differential  equations  are 
reduced  to  finite  difference  equations  and  solved  by  numeric 
methods . 
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Air  flow  through  a  test  cell  is  known  to  be  turbulent 
and  does  possess  regions  of  recirculation.   The  purpose  of 
this  study  was  to  apply  the  methods  described  by  Spalding 
to  the  flow  field  encompassing  the  engine  exhaust  and 
augmenter  tube  of  a  typical  Navy  jet  engine  test  cell. 

Flow  visualization  and  static  pressure  recovery  were 
primary  objectives  of  this  study.   Parametric  studies  were 
made  which  involved  the  variation  of  augmenter  geometry, 
engine  power  setting,  augmentation  ratio  and  cell  inlet 
flow  conditions. 

The  study  was  restricted  to  non-afterburning,  subsonic 
engine  exhaust  conditions. 
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II.   THEORETICAL  APPROACH 

Two-dimensional  fluid  flow  with  regions  of  reverse-flow 
can  be  described  mathematically  by  second-order  elliptic 
partial  differential  equations.   These  equations  can  be 
reduced  to  "finite-difference"  non-linear  algebraic  equations 
and  solved  simultaneously  by  an  iterative  procedure.   One 
such  procedure  is  detailed  by  Spalding  and  others  in  References 
1  and  2  and  was  the  method  used  in  this  study. 

The  conservation  of  mass,  momentum  and  energy  are  the 
fundamental  equations  used  in  the  problem  formulation.   The 
dependent  variables  chosen  by  Spalding  et  al.  [1,2]  are 
stream  functions  (ip)  and  vorticity  (w)  to  assure  mass  and 
momentum  conservation.   Temperature  was  used  as  the  dependent 
variable  for  energy  conservation  since  the  flow  was  considered 
inert.   Since  the  velocities  under  consideration  were  low, 
kinetic  heating  was  ignored.   Specific  heat  was  considered  to 
be  constant. 

Although  pressure  and,  to  a  large  extent,  velocity  are 
removed  from  the  formulation,  they  can  easily  be  "backed-out" 
of  the  vorticity  and  stream  functions  solutions  respectively. 

A  turbulence  model  is  needed  in  order  to  provide 
"effective"  transport  properties.   The  Jones-Launder 
turbulence  model  [2,3]  is  used  to  describe  the  variation 
of  the  effective  viscosity,  \i   ff.   Two  variables  are  used 
in  this  model:  turbulence  kinetic  energy,  K  and  turbulence 
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energy  dissipation  rate,  e.   The  model  is  a  time-averaged 
turbulence  model  and  the  two  associated  differential  equations 
can  be  cast  in  elliptic  form.   Thus,  the  solution  procedure 
for  tyt    w  and  T  also  applies  to  K  and  e .   The  relationship 
between  K  and  e  is: 


K2 

yeff  =  Cy  p  T 


where 


p   =   local  density 


C  =  empirically  determined  coefficient  =  0.09 


These  five  equations  together  with  the  equation  of  state 
for  a  perfect  gas  describe  the  complete  flow  field.   Choosing 
cylindrical  coordinates  (R,G,z)  allows  all  five  equations 
to  be  written  in  the  following  general  form: 

a      3(c.cf)) 


where  cj>  stands  for  \p ,    — ,  T,  K,  or  e.   The  functions  a,, 

b.,  c,  and  S,  are  given  in  Table  I. 

<p        <p  <p      •* 
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TABLE  I 
COEFFICIENT  VALUES 


|      R2     R2     yeff       s(o)/r) 


P  R 


1  -  H 

2  x  R 


yeff 
K        1      -|ii       1  -(G  -  pe) 

°K 
yeff  c^Ge    c2pe 


2 


a„  (   K   "   K   } 


yeff 
ah 


where 


F 


ef  f 
a„   =  Ti taken  to  be  1.0 


K,eff 


K,eff 


=  the  exchange  coefficient  of  the  kinetic 
energy  of  the  fluctuating  motion,  K. 
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eff 


e,ef  f 


taken  to  be  1.3 


the  exchange  coefficient  of  the  turbulence 
efeff     energy  dissipation  rate,  e. 


leff 


H  ~  r 


h,eff 


taken  to  be  1.0 


r.  cc  the  exchange  coefficient  of  heat 

h,ert 


=  V 


eff 


r  8V 


8V        8V        V   - 

2  +  <sir>2  +  <ir>2 


c,  and  c2  are  empirically  determined  coefficients. 


(c1  =  1.45,  c2  =  2.0) 


Neglecting  tangential  momentum,  the  vorticity  source  term 


S(w/R)   iS 


S(oJ/R)  -   "  R 


3z 


2      2 
V    +  V 

_z L.)  le. 

2     ;  8R 


2      2 

V    +  V 

3-  (I* 5_j  l£ 

8R  v     2     y  8z 


+  negligible  terms 
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The  first  term  of  the  general  equation  is  the  convective 
term  while  the  second  and  third  terms  are  the  diffusion  terms, 
The  last  term  is  the  source  term  which  includes  all  terms  not 
strictly  convective  or  diffusive. 

Once  the  equations  are  solved  and  a  converged  solution 
is  acquired,  the  static  pressure  distribution  in  the  flow 
field  can  be  calculated  by  integrating  the  momentum  equation 
in  the  axial  or  radial  direction.   The  resultant  equation  for 
the  gradient  in  the  axial  (z)  direction  is 


3P 

9z 


9_ 
9z 


eff 


r  9V 


9z 


div  V 


-  pvr 


'} 
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R3R  J  Ry 


eff 


9Vr 


+ 


9V 

9z" 


R 


-    RPVZVR 


} 


The  analogous  form  in  the  radial  (R)  direction  (neglecting 
6-momentum)  is 


1£   - 
9R 


where 


=  I  {-w.„(^ -§-,)] 


9z 


div  $ 


f      /  9Vz    8VR  \         1 
(^eff  (jtr+Jz-)-    PVzVR) 


dz 
9V 


s  h  {j%f f  [2  Xs  -  ! div  *]  -  rpvr2} 


8V 


z    +    1    9(RVR) 


9z  R 


9R 
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Pressure  is  then  obtained  by  the  following  simple 
integration  formula  (which  assumes  a  linear  variation  in 
pressure  gradient  between  grid  points) : 


,B  dP    ,    ~   1  r.3P.     ,3PX  ,  .  . 

/   87  dz   ~   2  t(a^>A  +  W  3  (ZB  "  ZA} 
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III.       MODEL    DISCUSSION 

A.  COORDINATE  SELECTION 

Although  the  test  cell  cross-section  is  basically  rectan- 
gular, the  engine  exit  and  augmenter  tube  can  be  accurately 
described  using  two-dimensional  axi-symmetrical  flow.   Of 
primary  interest  is  the  mixing  of  the  engine  exhaust  (primary 
flow)  and  ambient  air  (secondary  flow) .   This  process  pri- 
marily occurs  between  the  two  circular  "tubes"  and  within 
the  augmenter  tube.   Reference  4  suggests  that  test  cell 
wall  proximity  will  influence  primary/secondary  flow  mixing. 
However,  the  large  dimensions  of  the  cell  compared  with  the 
engine  and  augmenter  diameters  should  allow  the  use  of 
cylindrical  coordinates  with  little  error. 

B.  GEOMETRY 

Figure  1  is  a  profile  view  of  a  test  cell  located  at 
NARF  Alameda,  used  to  test  TF-41  engines.   Figure  2  is  the 
geometry  established  for  an  engine  test  while  Figure  3  is  the 
geometry  used  in  this  study.   The  actual  dimensions  are 
easily  varied  in  the  computer  program  allowing  parametric 
studies  to  be  made  to  assess  the  effects  of  geometry  change 
on  the  flow  field.   The  field  is  assumed  axi-symmetric  although 
in  practice  the  engine/augmenter  combination  is  located 
nearer  the  test  cell  floor  than  the  walls  or  ceiling.   The 
cell  ceiling  is  modeled  as  cylindrical  but  in  actuality  is 
quite  flat. 
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To  minimize  the  number  of  grid  points  required  in  the 
axial  (z)  direction,  the  cell  "inlet"  was  located  a  few  feet 
upwind  of  the  engine  exhaust. 

The  straight  augmenter  duct  "lip"  used  in  the  computer 
model  was  a  necessary  simplification  to  the  beveled  lip 
since  grid  points  in  the  coordinate  system  lie  along  mutually 
perpendicular  grid  lines.   Augmenter  tubes  in  Navy  test  cells 
are  by  no  means  uniform  in  design  and  inlet  shapes  vary 
greatly.   It  is  believed  that  a  straight  lip  design  does 
not  greatly  affect  the  flow  analysis. 

Augmenter  tubes  used  to  test  afterburner  operation 
usually  contain  quench  systems  and  flow  baffles.   These  were 
not  considered  in  this  analysis. 

C.   GRID  SPACING 

Various  grid  sizes  were  used  to  appraise  the  effect  of 
grid  size  on  the  predicted  flow  field.   The  largest  number 
of  grid  points  was  described  by  a  60x50  (zxR)  matrix  and  the 
smallest  by  a  20x25  matrix.   Regardless  of  the  size  of  the 
grid,  every  boundary  must  be  defined  by  a  grid  line.   In 
general,  grid  line  separation  should  be  smaller  as  one 
approaches  a  wall  since  gradients  in  properties  are  large 
in  this  region. 

Spalding  [1,2]  suggests  that  in  the  region  of  solid 
boundaries  each  succeeding  grid  spacing  not  differ  from  the 
last  by  more  than  a  factor  of  1.5. 
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D.   BOUNDARY  CONDITIONS 

All  dependent  variables  must  be  specified  in  some  manner 
on  the  boundaries.   The  boundaries  of  this  model  include: 
(1)   the  engine  exhaust  duct,  (2)  the  cell  inlet,  (3)  the 
solid  boundaries,  (4)  the  augmenter  outlet,  and  (5)  the 
geometric  centerline. 

1.   Engine  Exhaust  Duct 

a.  Stream  Function/Velocity 

Plug  flow  was  assumed  at  the  engine  exhaust 
plane.   The  axial  velocity  distribution  is  uniform  (Figure 
4)  to  produce  a  specified  mass  flow  and  no  radial  component 
of  velocity  exists. 

b.  Vorticity 

Since  no  velocity  gradients  exist,  vorticity 
is  identically  zero. 

c.  Mass  Fraction/Temperature 

A  non-dimensional  mass  fraction/temperature  was 
assigned  a  value  of  1.0  at  the  engine  exit.   This  corresponds 
to  the  actual  exhaust  temperature  existing  at  the  specified 
power  setting. 

d.  Turbulence  Kinetic  Energy 

Kinetic  energy  of  turbulent  motion  (K)  is  defined 
by 


K  -   I  (VR  +  V  +  V 


where   v,  ,  are  the  instantaneous  fluctuations 
in  free  stream  velocity. 
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For  isotropic  turbulence  the  expression  for  K 


3   2 
can  be  written  K  =  y  v  . 


Intensity  of  turbulence  is  defined  as 


1  ,  2  2  2, 

3  (vr  +  ve  +  vz  } 

V 


Assuming  isotropic  turbulence 

V 
where  V  is  the  mean  flow  velocity. 

Spalding  [1,2]  presents  an  empirical  relation 
for  K  in  uniform  flow: 


K  =  0.004  (V)2 


Setting 


3   2  —  2 

j   v    =   .004  (V)Z 


and  substituting  into  the  intensity  relation  we  find 


08  ^2 


I   =   ^ =   0.052 

V 
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Schlichting  [5]  experimentally  determined 
intensity  values  of  between  0.05  and  0.08  in  low  subsonic 
wind  tunnel  flow.   Thus,  it  appears  Spalding's  relation  is 
adequate  for  estimating  K  in  the  engine  exit  and  was  utilized 
as  the  boundary  condition. 

e.   Turbulence  Energy  Dissipation  Rate 

K  and  e  are  related  through  the  dissipation 
length  scale  I    [2]    by 


K3/2 

£        =   


In  this  program  t   was  taken  to  be: 

I     =   0.022  x  (Engine  Radius) 

2.   Inlet  to  Cell 

a.  Stream  Function/Velocity 

A  linear  variation  of  inlet  axial  velocity  was 
chosen  which  varied  from  a  positive  value  at  the  engine 
wall  to  zero  at  the  cell  ceiling  to  produce  a  specified 
secondary  mass  flow  (Figure  4) .   No  specific  information  was 
available  to  specify  this  flow  distribution  accurately. 
Lacking  more  specific  information,  the  radial  components  of 
velocity  (V_.)  were  fixed  at  zero. 

b.  Vorticity 

Vorticity  was  computed  using 

,   8V„    8V 
oj   _   !_  , R     _Zv 

R      R  l3z   "  8R  ' 
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c.  Mass-Fraction/Temperature 
Non-dimensional  temperature  and  mass-fraction 

were  assigned  a  value  of  0.0  at  the  cell  inlet.   This 
corresponded  to  an  antual  temperature  of  520  °R. 

d.  Turbulent  Kinetic  Energy  (K) 

K  was  specified  at  each  grid  point  using  the 
expression 


K   =   0.004  V2 


—2 

Where  V   is  the  local  flow  velocity 


e.   Turbulence  Dissipation 

C   =   K3/2/£ 

I      =   0.022  x  (Cell  Radius  -  Engine  Radius) 

3 .   Solid  Boundaries 

There  are  eight  solid  boundaries  described  in  the 
present  model.   Each  is  numbered  in  Figure  5  for  reference. 
Two  are  engine  surfaces,  two  are  cell  walls  and  four  describe 
the  augmenter  tube. 

a.   Velocity /Stream  Function 

Along  the  engine  surfaces  (1  and  2,  Figure  5) 
stream  function  is  constant  and  determined  from  the  specified 
mass  flow  through  the  engine.   The  remaining  solid  boundaries 
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have  the  stream  function  fixed  at  a  value  which  is  the  sum 
of  the  engine  mass  flow  and  cell  inlet  mass  flow.   All 
velocities  at  a  solid  boundary  are  zero, 
b.   Vorticity 

Vorticity  near  a  solid  boundary  was  computed  in 
one  of  two  ways  depending  on  the  magnitude  of  the  local  grid 
cell  Reynolds  number. 

(1)   Laminar  Sublayer  Case.   If  the  local  grid 
cell  Reynolds  number  (R  e  A^/yR)  is  found  to  be  less  than 
290  [1,2],  a  linear  velocity  gradient  is  assumed.   Using  the 
additional  assumption  that  the  velocity  gradient  in  the 
direction  perpendicular  to  the  wall  is  much  greater  than  the 
gradient  parallel  to  the  wall,  the  vorticity  is  calculated 
using  the  vorticity  definition  and  discarding  the  parallel 
term.   For  example,  given  a  wall  parallel  to  the  horizontal 
(z)  axis  with  flow  in  the  positive  z-direction,  the  vorticity 
would  be  found  by: 


03 


R     z  ~    _z_      _z  _   _p_ 

9z   '  3R       8R   '    AR       AR 

8VR 
where  ■* —  is  considered  small.   V  is  the  velocity  in  the 

O  Z  ir 

z-direction  at  node  P  (Figure  6) .   Since  velocity  at  the  wall 
is  zero,  it  is  easily  seen  that  V  is  twice  the  average 
velocity  between  the  wall  and  the  nearest  grid  point  off  the 
wall.   This  is  convenient  since  the  average  velocity  is 
related  to  the  stream  function  by 
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v 

V  P      1  Aif; 

zavg      2     PR  AR 


and  the  local  Reynolds  number  is 


Re 


y  R 
p 


Therefore,  the  vorticity  expression  is  given  by 

^z       ^P  2  Re  u 

»   -  -  AR      "  AR  j^f 

which  is  easily  calculated. 

(2)   Turbulent  Case.  When  the  local  grid  cell 

Reynolds  number  at  the  near  wall  node  exceeds  290  the  flow  is 

considered  turbulent.   Vorticity  is  then  defined  as 


Tshear 

co   =  

ywall 


where  j   ,  is  the  wall  shear  stress  and  u   ,,  is  the 

shear  Kwall 

viscosity  on  the  wall.   t  h    is  given  by  the  relation: 

Cf    -2 

Tshear  "  T  P  V 

Cf 
The  coefficient  (-—)    was  calculated  using 


Cf      .0284         .     _        A  if; 
-*-     =  5-f-  ,     where  Re   =  — -*- 

2       Re*2j  y  RP 
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(3)   Special  Points.   The  preceding  discussion 
of  vorticity  calculation  near  solid  boundaries  assumed  a 
well-developed  boundary  layer.   There  are  locations  in  the 
field  where  the  fluid  may  meet  a  parallel  solid  surface  and 
the  vorticity  calculation  becomes  suspect.   Such  a  point  is 
the  intersection  of  surface  5  and  surface  6  of  Figure  5 
(point  C) .   THere  are  two  problems  associated  with  this  point 
If  the  flow  is  predominantly  in  the  positive  z-direction,  a 
fully  developed  boundary  layer  will  not  exist  on  surface  5  at 
point  C.   If  the  flow  is  in  the  negative  R-direction  near 
surface  6,  the  identical  problem  arises.   Also,  each  surface 
(5  and  6)  will  have  a  unique  value  of  vorticity  at  point  C 
and  it  is  necessary  to  keep  them  properly  identified  in  the 
computer  program  (see  sketch  below) . 

•   •   • 

surface  5 


surface 
6 


R 
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No  completely  satisfactory  numerical  method 
exists  in  this  scheme  to  account  for  the  singularity.   Initial 
computer  runs  were  made  allowing  the  vorticity  at  point  C 
to  be  calculated  as  if  a  fully  developed  boundary  layer 
existed  on  both  surfaces  (5  and  6) .   Converged  solutions 
resulted  with  realistic  flow  visualizations  around  this  point 
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as  long  as  a  recirculation  region  did  not  develop  on  the 
augmenter  duct  lip  (surface  6)  . 

When  flow  developed  in  the  positive  R-direc- 
tion  near  surface  6  with  negative  z-direction  flow  on  surface 
5,  the  solution  would  not  converge  and  often  diverged  (see 
sketch  below) . 


•  *  0 

©   .; — -- 
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References  1  and  2  state  that  special 
attention  must  be  paid  to  points  which  define  enlargements 
in  pipe  flow  (Points  A  and  D  of  Figure  5)  .   In  order  to 
obtain  realistic  jet  spread,  stream  function  must  be  fixed 
on  the  node  immediately  downwind  of  these  points.   It  became 
evident  that  the  flow  would  usually  be  in  the  negative  z- 
direction  along  surface  5,  so  ty   was  fixed  at  the  node  circled 
in  the  sketch  above.   This  solved  the  divergence  problem 
occurring  under  certain  flow  conditions,  but  appeared  not  to 
alter  the  solution  appreciably.   Point  D  is  unusual  in  that 
it  defines  the  end  of  an  "infinitely  thin"  wall.   It  also 
defines  an  "enlargement"  and  for  this  reason,  the  immediate 
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downward  ip  value  was  fixed  at  the  value  assigned  to 
point  D. 


surface 
6 


surface  8 


surface 
7 


•  •  •   •  ip  fixed 
D 


dominant  flow 
direction 


For  purposes  of  calculating  the  vorticity  at  node  F  (see 
above  sketch) ,  the  vorticity  at  point  D  was  taken  to  be  the 
average  of  the  vorticities  computed  on  surface  6  and  7  shown 
in  the  above  sketch. 

c.   Turbulence  Kinetic  Energy 

Turbulent  flow  produces  steep  gradients  of  all 
transport  properties  near  walls.   To  reduce  the  required 
number  of  grid  points  and  to  reduce  computation  time  a 
"slip  value"  of  turbulence  kinetic  energy  (K)  is  calculated 
on  the  wall.   The  "slip  value"  is  fixed  to  produce  an  arith- 
fnetic-mean  gradient  between  the  wall  and  adjacent  node,  which 
equals  the  gradient  of  the  Couette  flow  solution  midway 
between  the  nodes.   A  complete  analysis  of  this  method  can 
be  found  in  References  1  and  6. 
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d.   Turbulence  Dissipation 

The  effective  boundary  for  turbulence  dissipa- 
tion (e)  is  taken  to  be  the  n  grid  points  adjacent  to  the 
wall.   Wolfshtein  [6]  presents  the  value  at  the  near-wall 
node  as : 


E 


=   (!&)  (Kp)3/2/AR 


where   e    =  near-wall  value 


c   =  empirical  constant  =   0.09 

a   =   empirical  constant   =   0.22 

Kp   =   near-wall  turbulence  kinetic 
energy  value 

AR  =  wall  to  near-wall  node 
distance 


e.   Temperature/Mass  Fraction 

All  solid  boundaries  were  assumed  adiabatic 
and,  of  course,  impermeable.   Therefore,  no  gradients  in 
temperature  or  mass  fraction  existed  between  the  walls  and 
near-wall  nodes. 

4 .  Augmenter  Exit 

Axial  gradients  of  all  dependent  variables  were 
taken  as  zero. 

5.  Geometric  Centerline 

Radial  gradients  of  w/R,  K,  c,  T  and  mass  fraction 
are  taken  as  zero.   Also  the  centerline  is  the  stream  function 
(ijj  =  0)  reference  line. 
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E.  INITIAL  CONDITIONS 

The  flow  field  must  have  initial  values  assigned  to  all 
variables.   Experience  is  required  to  fill  the  field  with 
compatible  variables  to  prevent  divergence  in  the  first  few 
iterations . 

It  was  found  that  a  reasonably  accurate  stream  function 
distribution  was  essential  to  prevent  divergence.   In  par- 
ticular, the  predominant  flow  direction  must  be  established 
and  large  discontinuities  in  stream  function  should  not 
exist.   The  flows  in  the  augmenter  tube  and  the  free  stream 
area  between  the  engine  exhaust  and  augmenter  tube  were 
initialized  as  well-developed  pipe  flow.   This  allowed  com- 
patible vorticity  values  to  be  calculated  in  these  regions. 

Density,  mass  fraction,  viscosity,  and  turbulence 
kinetic  energy  were  specified  using  reference  values  taken 
from  either  the  cell  inlet  or  engine  exit.   Turbulence 
dissipation  was  calculated  using  the  relationship  discussed 
under  boundary  conditions. 

F.  SOLUTION  PROCEDURES 

To  counter  computational  instabilities  in  the  first  400 
iterations,  stream  function  variation  between  iterations  was 
under-relaxed  by  a  factor  of  0.3.   All  other  variables  were 
under-relaxed  by  0.5.   The  relaxation  parameters  were 
increased  to  0 . 6  and  0.8,  respectively,  after  400  iterations 

•  Various  flows,  geometries,  and  grids  were  investigated 
successfully  using  the  above  relaxation  parameters.   They 
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precluded  divergence  in  the  most  unstable  flow  condition 
handled:   high  engine  mass  flow  with  low  augmentation  ratio 
(200  lbm/sec,  0.25,  respectively).   However,  it  should  be 
noted  that  initial  conditions  must  be  well  established  for 
these  relaxation  parameters  to  be  effective. 
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IV.   MODEL  DEVELOPMENT 

The  geometry  shown  in  Figure  3a  was  chosen  to  represent 
a  test  cell  section  containing  an  engine  exhaust  nozzle  and 
augmenter  tube.   Approximate  dimensions  were  taken  from  a 
TF-41  turbofan  engine  undergoing  testing  at  NARF  Alameda. 

A.   GRID  SPACING 

A  grid  composed  of  twenty  lines  in  the  z-direction  and 
twenty-five  lines  in  the  R-direction  (20x25)  was  chosen  to 
initially  describe  the  flow  field.   (Intersections  of 
mutually  perpendicular  grid  lines  define  grid  points . ) 
This  number  of  grid  lines  allowed  the  use  of  small  grid 
spacing  near  solid  walls  but  produced  relatively  large 
spacing  away  from  walls. 

Computation  time  is  directly  proportional  to  the  number 
of  grid  points  in  the  field.   A  small  number  of  grid  lines 
was  chosen  to  minimize  computer  time  during  program  debugging 
Three  hundred  iterations  typically  required  7.5  minutes 
(excluding  compile  time) .   It  should  be  noted  that  vorticity 
source  terms  were  ignored  during  this  initial  period. 

Vorticity  source  terms  and  pressure  recovery  calculations 
were  added  only  after  the  program  was  successfully  iterating 
to  a  converged  solution.   The  number  of  grid  lines  was 
increased  to  40x40  (zxR)  to  allow  for  pressure  gradient 
calculations  without  the  use  of  wall  or  near-wall  variable 
values.   Otherwise,  large  errors  accumulated  in  the  pressure 
integration  process. 
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Eleven  hundred  iterations  typically  required  109  minutes 
of  computation  time  (excluding  compile  time) .   This  time 
includes  one  pressure  calculation  along  eight  lines  in  the 
flow  field  and  those  computations  performed  in  the  flow 
initialization  process.   However,  the  total  time  for  the 
latter  calculations  amounts  to  less  than  one  minute. 

B.  CONVERGENCE    CRITERIA 

Convergence  rate  was  influenced  by  flow  conditions,  grid 
spacing  and  relaxation  parameters.   Two  criteria  were  used 
to  appraise  the  convergence  level  of  the  solution. 

First,  the  percentage  change  in  u,  ty,    K,  e  and  mass 
fraction  (at  all  grid  points  in  the  field)  must  be  arbi- 
trarily small  between  successive  iterations.   This  was  done 
computationally  by  storing  the  largest  change  occurring  in 
the  field  for  each  variable  and  printing  that  value  after 
each  iteration. 

The  second  convergence  criterion  proved  to  be  the  most 
conclusive.   That  is,  when  pressure  calculations  along  each 
grid  line  remained  unchanged  (with  increased  iterations) 
convergence  was  assumed. 

C.  PLOTTING  PACKAGE 

Solution  values  of  stream  function  and  temperature  were 
transferred  to  computer  cards  at  the  conclusion  of  each  run. 
These  cards  were  then  used  as  input  to  a  second  program 
which  plotted  streamline  contours  and  lines  of  constant 
temperature  on  separate  graphs.   A  CALCOMP  Model  765  plotter 
was  utilized  for  this  purpose. 
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V.   PARAMETRIC  STUDY 

Geometry  and  inlet  conditions  were  varied  to  appraise 
model  utility  and  sensitivity  to  flow  conditions.   TF-41 
engine  test  data  were  used  to  establish  engine  exhaust  static 

temperature  and  velocity.   Static  pressure  at  the  engine 

2 
exit  was  fixed  at  atmospheric  pressure  (2116.0  lbf/ft  )  for 

this  study.   Engine  exit  velocity  was  found  by  computing  p 

from  the  equation  of  state  for  a  perfect  gas  and  using  the 

one-dimensional  mass  conservation  equation. 


P 

RT 


and 


V 


m 
engine 

engine    p  A 

3      K   engine 


where  m  =   total  mass  flow  through  engine 
(including  bypass  air) 

A  =   area  of  exhaust  plane 


Engine  exit  static  temperature  was  computed  using  the 
turbine  outlet  temperature  and  the  bypass  ratio  given  in 
Reference  7.   A  summary  of  these  parameters  is  shown  in 
Table  II. 
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TABLE  II 
TF-41  TEST  CONDITIONS 


Setting  *   (lb™/sec>  Texit    (°R»        Vexit(f ^sec) 


Exit 
Mach 
Number 


Idle  100  800  594  0.43 

80%  150  900 

90%  200  927         1376  0.92 

94%  238  955 

97%  249  970 

100%  263  1001 


A  specific  flow  condition  was  chosen  to  provide  a  solution 
to  which  all  succeeding  solutions  could  be  compared.   That 
condition  was:   idle  engine  power  setting  with  an  augmentation 
ratio  (AR)  of  0.5.   The  geometry  was  fixed  as  shown  in 
Figure  3. 

The  "reference"  flow  condition  was  chosen  as  it  provided 
low  subsonic  flow  conditions  at  the  engine  exhaust  plane 
(M  =  0.4  3)  and  the  augmentation  ratio  was  typical  for  actual 
test  conditions. 
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Pressure  calculations  were  made  along  the  dashed  lines 
shown  in  Figure  7.   The  lines  show  only  approximate  locations 
since  grid  spacing  was  often  altered  when  geometry  was 
varied. 

Eight  computer  runs  were  made  to  examine  the  usefulness 
of  the  model.   Each  run  was  allowed  approximately  2.5  hours 
(1400  iterations)  to  assure  solution  convergence  although 
most  runs  achieved  practical  convergence  prior  to  this 
time.   A  summary  of  the  study  follows  in  Table  III. 

TABLE  III 
FLOW  CONDITIONS 

1.  Reference  flow  condition  was  established 
as  discussed  above:   Idle  Power,  AR  =  0.5, 
Figure  3a  geometry. 

2.  Two  runs  were  made  varying  the  AR:   Idle 
Power,  AR  =  0.75  and  1.0,  Figure  3a 
geometry. 

3.  Two  runs  were  made  to  change  the  length 
of  the  augmenter  duct  lip:  Idle  Power, 
AR  =  0.5,  Figure  3b,  geometries. 

4 .  One  run  varied  the  engine  to  augmenter 
separation  distance:   Idle  Power,  AR  =  0.5, 
Figure  3c,  geometry. 

5.  One  run  varied  the  augmenter  radius: 

Idle  Power,  AR  =  0.5,  Figure  3d,  geometry. 

6.  One  run  changed  the  power  setting  of  the 
engine:   90%  power,  AR  =  1.0,  Figure  3a 
geometry. 
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VI.   RESULTS  AND  DISCUSSION 

A.   FLOW  CONDITION  1  (Reference  Condition) 

Streamlines  and  static  temperature  contours  are  shown 
in  Figures  8a  and  8b,  respectively.   Pressure  gradients  along 
axial  and  radial  grid  lines  (Figure  7)  are  graphically 
displayed  in  Figures  8c  through  8f.   Pressure  contours  along 
axial  and  radial  grid  lines  are  displayed  in  Figures  8g 
through  8i.   Figure  8j  presents  the  variation  in  augmenter 
pressure  rise  as  a  function  of  iteration  number. 

1 .   Streamlines/Temperature  Contours 

Two-dimensional  flow  solutions  are  produced  by  this 
study.   Therefore,  the  flow  field  consists  of  a  one-radian 
"portion"  of  the  axisymmetric  field  (Figure  8a)  and  the 
streamlines  are  labeled  accordingly.   For  example,  the  mass 
flow  rate  of  a  TF-41  engine  is  100  lb  /sec  at  idle  power 
setting.   The  two-dimensional  mass  flow  would  be 


actual 
2-D      2tt 


100  lb  /sec 
m 


=   15.9  lb  /sec-radian 


2  IT  ITV 

as  is  shown  in  Figure  8a. 

Two  distinct  regions  of  recirculation  are  apparent 
in  Figure  8a.   One  prominent  region  exists  above  the  augmenter 
tube  bounded  on  two  sides  by  the  test  cell  ceiling  and  wall. 
Points  of  stagnation  are  shown  located  on  the  exterior 
augmenter  tube  surface  and  cell  wall. 
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The  second  prominent  region  of  recirculation  is 
located  downwind  of  the  augmenter  inlet  lip.   Close  exam- 
ination of  the  streamline  having  its  origin  at  the  lip 
reveals  axial  flow  for  a  short  distance  followed  by  rapid 
directional  change.   Axial  flow  could  have  been  the  result 
of  fixing  the  stream  function  at  the  first  downwind  corner 
as  discussed  in  Section  III .D. 3 .b . ( 3) .   However,  an  additional 
computer  run  was  made  with  the  same  flow  conditions  except 
that  the  stream  function  at  the  downwind  node  was  allowed  to 
be  computed  normally.   The  solution  for  stream  function  was 
essentially  identical  to  the  first  run.   In  any  case,  the 
lip  has  only  small  effects  on  the  shape  and  size  of  the 
recirculation  region.   The  recirculation  region  is  formed 
at  the  low  thrust,  low  augmentation  ratio  condition  pri- 
marily by  the  radial  component  of  the  augmentation  air  as 
it  flows  into  the  augmenter. 

The  region  of  flow  reattachment  is  well  defined  on 
the  augmenter  interior  surface  and  is  located  approximately 
two  augmenter  diameters  downstream  of  the  inlet. 

Temperature  contours  shown  in  Figure  8b  are  an 
indication  of  the  amount  of  mixing  occurring  between 
engine  flow  and  secondary  air.   From  the  non-dimensional 
mass  fraction  and  temperature  definitions,  a  thoroughly 
mixed  flow  would  have  a  static  temperature  of  707  °R. 
Under  these  flow  conditions,  a  well-mixed  flow  occurs 
approximately  four  augmenter  diameters  downstream 
of  the  augmenter  inlet.   The  temperature  contours 


42 


clearly  show  the  effect  of  secondary  air  entrainment  and 
recirculation  within  the  augmenter  tube  on  the  temperature 
field. 

2  .   Pressure  Recovery 

a.   Axial  (z-direction) 

Pressure  gradients  along  lines  A-A  through  E-E 
(Figure  7)  are  shown  in  Figures  8c,  d  and  e.  Lines  labeled 
A-A  and  B-B  are  axial  lines  having  their  origin  at  the 
engine  exit  plane  and  terminating  one  node  upstream  of  the 
augmenter  exit.  Lines  C-C,  D-D  and  E-E  extend  between  the 
third  node  in  the  axial  direction  and  the  node  immediately 
upstream  of  the  augmenter  exit. 

Large  negative  pressure  gradients  exist  along  all 
lines  from  the  engine  exit  plane  to  just  inside  the  augmenter 
inlet  plane.   Lines  A-A  and  B-B  exhibit  the  largest  gradients 
in  this  region  and  examination  of  the  solution  reveals  the 
exhaust  jet  rapidly  departing  from  the  uniform  inlet  velocity 
(plug  flow)  condition.   Secondary  air  is  being  entrained  in 
this  area  with  line  B-B  being  nearest  the  turbulent  mixing 
zone  of  the  exhaust  jet.   This  may  account  for  the  largest 
pressure  gradients  existing  along  line  B-B. 

Pressure  gradients  along  axial  line  F-F  (Figure  7) 
are  very  small  and  are  not  shown.   Static  pressures  along 
axial  lines  A-A  through  E-E  are  shown  in  Figures  8h  and  i. 
These  pressure  profiles  are  all  initialized  at  2116.0 
lbf/ft2. 
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b.  Radial  (R-direction) 

Pressure  gradients  and  static  pressures  along 
G-F  are  shown  in  Figures  8f  and  8g,  respectively.   The  engine 
exhaust  jet  turbulent  mixing  zone  is  clearly  indicated  by 
fluctuating  pressure  gradients  (Figure  8f) .   Static  pressure 
is  seen  to  increase  toward  the  engine  exhaust  jet  centerline 
which  is  due  primarily  to  the  secondary  air  momentum  in 
the  negative  R-direction. 

No  appreciable  pressure  gradients  existed  along 
radial  line  H-H  which  is  located  downstream  of  the  flow 
reattachment  point. 

c.  Total  Axial  Static  Pressure  Rise 

Total  rise  (AP)  was  defined  as  the  static  pressure 
difference  between  end-points  of  the  lines  discussed  above. 
AP  is  plotted  against  number  of  iterations  in  Figure  8 j .   The 
upwind  end  pressure  for  all  lines  is  reference  pressure  but 
each  line  is  seen  to  exhibit  a  unique  value  for  AP  at  1400 
iterations.   No  pressure  gradients  exist  in  the  radial 
direction  downwind  of  the  reattachment  point.   Therefore, 
AP  values  should  be  identical  in  magnitude  for  lines  A-A 
through  E-E  when  the  solution  is  converged. 

Close  examination  of  the  pressure  gradients 
(Figures  8c  through  8e)  reveals  line  E-E  to  have  small,  "well 
behaved"  gradients  compared  with  lines  A-A  through  D-D. 
Figure  8j  indicates  that  the  AP  calculation  for  line  E-E 
"converged"  relatively  early  (1000  iterations) .   For  these 
reasons,  pressure  rise  was  taken  to  be  that  computed  for 
E-E  (47  lbf/ft2) . 
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The  differences  in  the  pressure  rise  calculated 
along  axial  grid  lines  can  be  due  to  several  things.   The 
number  of  iterations  has  a  dominant  effect  on  the  calculated 
pressure  rise.   The  dependent  variables  (co,  \p ,    K,  e ,  T)  vary 
only  slightly  with  increased  iteration  time.   However,  the 
pressure  calculation  is  apparently  very  sensitive  to  small 
changes  in  the  stream  function  and  density.   This  is  probably 
because  —■   is  a  function  of  second  derivatives  of  ty.      In 

dZ 

addition,  pressure  rise  is  calculated  by  assuming  linear 
variations  in  gradient  between  grid  points  and  pressure  has 
been  held  fixed  for  purposes  of  calculating  the  density 
(i.e.,  p  =  p(T)).   Improved  accuracy  might  be  obtained  if 
pressure  were  calculated  periodically  in  the  solution  in 
order  that  pressure  effects  on  density  may  be  taken  into 
account.   It  was  found  that  finer  grid  spacing  did  not 
appreciably  affect  the  pressure  calculation.   Another  possible 
error  source  may  enter  through  the  use  of  upwind  differences 
and/or  the  inadequacy  of  the  Jones-Launder  model  for  the  flow 
geometry  under  investigation.   However,  checks  on  the 
solution  indicated  that  mass  was  conserved  at  each  axial 
location  with  less  than  14  percent  error. 

B.   FLOW  CONDITION  2 

Augmentation  ratio  (AR)  was  varied  with  idle  power  setting 
to  assess  the  effect  on  total  static  pressure  rise  through  the 
system.   Solution  convergence  and  pressure  gradients  were 
similar  to  Flow  Condition  1  and  axial  line  E-E  was  again 
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chosen  to  determine  the  value  of  AP.   The  solutions  along 
the  other  axial  lines  had  approximately  the  same  slope  on  an 
AR  vs.  AP  plot  but  were  different  in  magnitude  as  discussed 
above. 

Figure  9a  shows  decreasing  AP  with  increased  AR.   This 
result  qualitatively  correlates  with  Bailey  [8,9]  who  experi- 
mentally varied  exit  pressure  on  a  small  jet  pump  and  deter- 
mined secondary  mass  entrainment.   As  the  augmentation  air 
is  increased  for  a  fixed  engine  flow  there  is  less  momentum 
differential  between  the  two  streams.   This  apparently  causes 
the  decreased  pressure  rise  with  increased  AR.   In  actual 
test  cell  operation  there  is  a  unique  AR  for  each  engine 
flow  rate  (assuming  the  augmenter  is  fixed  in  design  and 
position) .   The  actual  value  is  determined  from  the  pressure 
rise  in  the  augmenter  tube  being  equal  to  the  pressure  losses 
through  the  rest  of  the  test  cell.   If  pressure  losses 
through  the  inlet  and  exit  stacks  (with  their  associated 
acoustic  and  pollution  abatement  equipment)  and  through  the 
horizontal  portion  of  the  cell  can  be  estimated  as  a  function 
of  flow  velocity,  then  the  AR  can  be  found  in  an  iterative 
manner  by  using  the  AP  vs.  AR  plot  determined  with  the  model 
for  a  particular  augmenter  design  and  position. 

Figures  9b  and  9c  show  the  effect  of  AR  on  the  stream 
function  and  temperature  distributions,  respectively.   As 
the  AR  is  increased  the  axial  velocity  component  of  the 
secondary  air  dominates  over  the  radial  component.   The  latter 
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component  was  observed  above  to  be  the  major  cause  of  the 
recirculation  region.   Thus,  as  the  AR  is  increased,  the 
recirculation  region  within  the  augmenter  is  seen  to  decrease 
in  size  and  mixing  occurs  more  rapidly. 

C.  FLOW  CONDITION  3 

The  augmenter  duct  lip  length  was  varied  by  0.6  3  ft 
(2.1  ft  to  2.73  ft  augmenter  inlet  radius).   Stream  function 
and  temperature  fields  are  shown  in  Figures  10a  through  lOd. 
No  appreciable  change  in  calculated  pressure  rise  was  observed 
This  result  apparently  occurs  because  at  low  augmentation 
ratios  the  size  of  the  inlet  lip  has  very  little  effect  on 
the  size  of  the  recirculation  region  within  the  augmenter. 
A  greater  effect  of  duct  lip  length  on  pressure  rise  could 
be  expected  at  higher  augmentation  ratios . 

D.  FLOW  CONDITION  4 

Engine  to  augmenter  inlet  separation  was  increased  to 
2.5  ft  (from  1.5  ft).   Stream  function  and  temperature 
contours  are  included  as  Figures  11a  and  lib.   No  appreciable 
change  in  pressure  rise  was  noted  along  axial  grid  line  E-E 
when  compared  with  Flow  Condition  1. 

The  stream  function  solution  indicated  that  no  appreciable 
change  occurred  within  the  augmenter  when  the  distance  was 
increased.   Again,  a  more  significant  effect  could  be  expected 
for  higher  augmentation  ratios.   Moving  the  augmenter  closer 
to  the  engine  should  have  a  major  effect  on  the  recirculation 


47 


region  and  the  calculated  pressure  rise.   As  discussed  above, 
the  inlet  velocity  profile  for  the  augmentation  air  was 
assumed  to  increase  linearly  from  the  cell  ceiling  to  the 
engine  at  a  location  just  upstream  of  the  engine  exhaust. 
This  profile  would  change  as  the  augmenter  position  is  changed 
Since  test  cells  are  designed  to  provide  uniform  flow  at  the 
engine  inlet,  future  investigations  with  the  model  should  be 
made  with  the  secondary  air  inlet  specified  as  uniform  at  the 
engine  inlet  plane.   In  this  case,  the  velocity  profile  near 
the  engine  exhaust  plane  would  be  free  to  develop  in  a 
realistic  manner. 

E.  FLOW  CONDITION  5 

Augmenter  radius  was  increased  by  0.5  ft  to  3.5  ft  while 
keeping  the  duct  inlet  diameter  fixed.   Figures  12a  and 
12b  present  the  temperature  and  stream  function  contours. 

Pressure  rise  (AP)  along  E-E  was  diminished  to  37  lbf/ft2 

2 
from  the  reference  condition  value  of  47  lbf/ft  .   Although  ' 

only  one  run  was  made,  this  pressure  loss  is  considerable  an 

indicates  AP  decreases  with  increasing  augmenter  diameter 

(given  fixed  inlet  flow  conditions) .   The  increased  size  of 

the  recirculation  region  (which  consists  primarily  of 

augmentation  air)  may  contribute  to  this  pressure  loss. 

F.  FLOW  CONDITION  6 

Engine  power  setting  was  increased  to  90%  and  the  A  R 
was  set  at  1.0  which  is  typical  for  an  actual  engine  test. 
Figure  13a  shows  the  recirculation  region  smaller  with  flow 


48 


reattachment  farther  upstream.   The  temperature  contours 
(Figure  13b)  show  engine  exhaust  and  secondary  air  being 
thoroughly  mixed  farther  down  the  augmenter  tube  than 
occurred  at  idle  flow. 

Pressure  rise  versus  number  of  iterations  is  shown  in 
Figure  13c.   The  numerical  values  of  the  dependent  variables 
did  not  change  appreciably  after  1400  iterations.   However, 
AP ' s  differed  in  value  considerably  along  lines  A-A  through 
E-E,  so  no  definitive  pressure  rise  could  be  derived  from  the 
results.   The  solution  procedure  (which  assumes  incompressible 
flow  and  uses  the  associated  elliptic  equations)  used  in 
this  study  is  questionable  when  applied  at  the  higher  subsonic 
Mach  numbers . 
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/ 
VII.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  model  provides  two-dimensional  flow  visualizations 
in  turbojet  test  cells  and  augmenter  tubes  for  low  subsonic 
flow  conditions.  These  can  be  used  to  identify  regions  of 
recirculation  and  to  assess  the  amount  of  mixing  occurring 
between  engine  exhaust  gases  and  secondary  air.  Optimum 
locations  for  pollution  sampling  equipment  can  be  selected 
by  examining  the  numerical  solutions. 

Static  pressure  rise  within  the  augmenter  tube  can  be 
calculated  for  a  given  flow  condition  and  augmenter  design 
from  converged  solutions,  provided  that  regions  of  "steep" 
gradients  in  variables  are  avoided.   Total  pressure  rise 
through  the  system  can  be  found  and  factors  which  influence 
this  rise  (augmenter  design  and  location,  engine  flow  rate, 
test  cell  size,  and  etc.)  can  be  studied. 

High  subsonic  flows  invalidate  the  use  of  partial 
differential  equations  of  the  elliptic  type;  therefore, 
the  model  provides  realistic  solutions  only  at  low  subsonic 
Mach  numbers. 

B.  RECOMMENDATIONS 

1.   Additional  computer  runs  should  be  made  to  firmly 
establish  the  effects  of  augmenter  and  cell  geometry  on  flow 
mixing  and  pressure  rise  (AP) .   For  example,  the  effects  of 
moving  the  augmenter  inlet  closer  to  the  engine  exit  plane 
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and  moving  the  test  cell  wall  forward  to  the  augmenter  inlet 
plane  should  be  considered.   In  addition,  the  program  could 
also  be  easily  modified  to  extend  the  augmenter  duct  lip  in 
the  opposite  direction. 

2.  The  cell  "inlet"  position  should  be  moved  upstream 
to  the  engine  inlet  plane  by  the  addition  of  a  few  grid 
lines  in  the  z-direction.   Plug  flow  could  be  assumed  at 
the  "new"  cell  inlet  and  the  resultant  flow  field  would  show 
the  actual  velocity  distribution  at  the  "cell  inlet"  used 

in  this  study.   In  addition,  this  would  permit  calculation 
of  exhaust  gas  recirculation. 

3.  Static  pressure  and  velocity  measurements  should  be 
taken  under  various  test  conditions  at  NARF  Alameda.   The 
size  and  shape  of  the  recirculation  region  located  above 
the  augmenter  should  be  qualitatively  verified  by  photo- 
graphing a  tufted  grid  placed  between  the  test  cell  augmenter 
tube  and  ceiling.   These  data  would  help  in  validation/ 
improvement  of  the  computer  model. 
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FIGURE  3b.   COMPUTER  MODEL  GEOMETRY:   FLOW  CONDITION  3 
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FIGURE  3c.   COMPUTER  MODEL  GEOMETRY:   FLOW  CONDITION  4 
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FIGURE  3d.   COMPUTER  MODEL  GEOMETRY:  FLOW  CONDITION  5 
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FIGURE  4.     VELOCITY  PROFILES   IN  INLET  FLOWS 
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FIGURE  5.  IDENTIFICATION  OF  BOUNDARIES 
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FIGURE  8f.      RADIAL  PRESSURE  GRADIENTS 
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